[INFO   ] Starting gmx_MMPBSA v1.5.2
[DEBUG  ] WDIR          : /media/romi/My Book Duo/completed-simulation-oracle/auto-immune/rs2/5mv4-rs2/312k/free-energy
[DEBUG  ] AMBERHOME     : /home/romi/anaconda3/envs/gmxMMPBSA
[DEBUG  ] PYTHON EXE    : /home/romi/anaconda3/envs/gmxMMPBSA/bin/python
[DEBUG  ] PYTHON VERSION: 3.9.12 (main, Apr  5 2022, 06:56:58) [GCC 7.5.0]
[DEBUG  ] MPI           : /home/romi/anaconda3/envs/gmxMMPBSA/bin/mpirun
[DEBUG  ] ParmEd        : 3.4.3
[DEBUG  ] OS PLATFORM   : Linux-5.15.0-58-generic-x86_64-with-glibc2.31
[DEBUG  ] OS SYSTEM     : Linux
[DEBUG  ] OS VERSION    : #64~20.04.1-Ubuntu SMP Fri Jan 6 16:42:31 UTC 2023
[DEBUG  ] OS PROCESSOR  : x86_64

[INFO   ] Command-line
  gmx_MMPBSA -O -i mmpbsa.in -cs md_0_10.tpr -ci index.ndx -cg 18 19 -ct dynamic-nopbc.xtc -cp topol.top -o FINAL_RESULTS_MMPBSA.dat -eo FINAL_RESULTS_MMPBSA.csv

[DEBUG  ] |Input file:
[DEBUG  ] |--------------------------------------------------------------
[DEBUG  ] |Input file generated by gmx_MMPBSA (v1.5.2)
[DEBUG  ] |Be careful with the variables you modify, some can have severe consequences on the results you obtain.
[DEBUG  ] |
[DEBUG  ] |# General namelist variables
[DEBUG  ] |&general
[DEBUG  ] |  sys_name             = ""                                             # System name
[DEBUG  ] |  startframe           = 10                                             # First frame to analyze
[DEBUG  ] |  endframe             = 9900                                       # Last frame to analyze
[DEBUG  ] |  interval             = 1600                                         # Number of frames between adjacent frames analyzed
[DEBUG  ] |  forcefields          = "oldff/leaprc.ff99SB,leaprc.gaff"              # Define the force field to build the Amber topology
[DEBUG  ] |  ions_parameters      = 1                                              # Define ions parameters to build the Amber topology
[DEBUG  ] |  PBRadii              = 3                                              # Define PBRadii to build amber topology from GROMACS files
[DEBUG  ] |  temperature          = 312                                        # Temperature
[DEBUG  ] |  qh_entropy           = 0                                              # Do quasi-harmonic calculation
[DEBUG  ] |  interaction_entropy  = 0                                              # Do Interaction Entropy calculation
[DEBUG  ] |  ie_segment           = 25                                             # Trajectory segment to calculate interaction entropy
[DEBUG  ] |  c2_entropy           = 0                                              # Do C2 Entropy calculation
[DEBUG  ] |  assign_chainID       = 0                                              # Assign chains ID
[DEBUG  ] |  exp_ki               = 0.0                                            # Experimental Ki in nM
[DEBUG  ] |  full_traj            = 0                                              # Print a full traj. AND the thread trajectories
[DEBUG  ] |  gmx_path             = ""                                             # Force to use this path to get GROMACS executable
[DEBUG  ] |  keep_files           = 2                                              # How many files to keep after successful completion
[DEBUG  ] |  netcdf               = 0                                              # Use NetCDF intermediate trajectories
[DEBUG  ] |  solvated_trajectory  = 1                                              # Define if it is necessary to cleanup the trajectories
[DEBUG  ] |  verbose              = 1                                              # How many energy terms to print in the final output
[DEBUG  ] |/
[DEBUG  ] |
[DEBUG  ] |# (AMBER) Generalized-Born namelist variables
[DEBUG  ] |&gb
[DEBUG  ] |  igb                  = 5                                              # GB model to use
[DEBUG  ] |  intdiel              = 1.0                                            # Internal dielectric constant for sander
[DEBUG  ] |  extdiel              = 78.5                                           # External dielectric constant for sander
[DEBUG  ] |  saltcon              = 0.0                                            # Salt concentration (M)
[DEBUG  ] |  surften              = 0.0072                                         # Surface tension
[DEBUG  ] |  surfoff              = 0.0                                            # Surface tension offset
[DEBUG  ] |  molsurf              = 0                                              # Use Connelly surface ('molsurf' program)
[DEBUG  ] |  msoffset             = 0.0                                            # Offset for molsurf calculation
[DEBUG  ] |  probe                = 1.4                                            # Solvent probe radius for surface area calc
[DEBUG  ] |  ifqnt                = 0                                              # Use QM on part of the system
[DEBUG  ] |  qm_theory            = ""                                             # Semi-empirical QM theory to use
[DEBUG  ] |  qm_residues          = ""                                             # Residues to treat with QM
[DEBUG  ] |  qmcharge_com         = 0                                              # Charge of QM region in complex
[DEBUG  ] |  qmcharge_lig         = 0                                              # Charge of QM region in ligand
[DEBUG  ] |  qmcharge_rec         = 0                                              # Charge of QM region in receptor
[DEBUG  ] |  qmcut                = 9999.0                                         # Cutoff in the QM region
[DEBUG  ] |  scfconv              = 1e-08                                          # Convergence criteria for the SCF calculation, in kcal/mol
[DEBUG  ] |  peptide_corr         = 0                                              # Apply MM correction to peptide linkages
[DEBUG  ] |  writepdb             = 1                                              # Write a PDB file of the selected QM region
[DEBUG  ] |  verbosity            = 0                                              # Controls the verbosity of QM/MM related output
[DEBUG  ] |  alpb                 = 0                                              # Use Analytical Linearized Poisson-Boltzmann (ALPB)
[DEBUG  ] |  arad_method          = 1                                              # Selected method to estimate the effective electrostatic size
[DEBUG  ] |/
[DEBUG  ] |
[DEBUG  ] |# (AMBER) Possion-Boltzmann namelist variables
[DEBUG  ] |&pb
[DEBUG  ] |  ipb                  = 2                                              # Dielectric model for PB
[DEBUG  ] |  inp                  = 2                                              # Nonpolar solvation method
[DEBUG  ] |  sander_apbs          = 0                                              # Use sander.APBS?
[DEBUG  ] |  indi                 = 1.0                                            # Internal dielectric constant
[DEBUG  ] |  exdi                 = 80.0                                           # External dielectric constant
[DEBUG  ] |  emem                 = 4.0                                            # Membrane dielectric constant
[DEBUG  ] |  smoothopt            = 1                                              # Set up dielectric values for finite-difference grid edges that are located across the solute/solvent dielectric boundary
[DEBUG  ] |  istrng               = 0.0                                            # Ionic strength (M)
[DEBUG  ] |  radiopt              = 1                                              # Use optimized radii?
[DEBUG  ] |  prbrad               = 1.4                                            # Probe radius
[DEBUG  ] |  iprob                = 2.0                                            # Mobile ion probe radius (Angstroms) for ion accessible surface used to define the Stern layer
[DEBUG  ] |  sasopt               = 0                                              # Molecular surface in PB implict model
[DEBUG  ] |  arcres               = 0.25                                           # The resolution (Å) to compute solvent accessible arcs
[DEBUG  ] |  memopt               = 0                                              # Use PB optimization for membrane
[DEBUG  ] |  mprob                = 2.7                                            # Membrane probe radius in Å
[DEBUG  ] |  mthick               = 40.0                                           # Membrane thickness
[DEBUG  ] |  mctrdz               = 0.0                                            # Distance to offset membrane in Z direction
[DEBUG  ] |  poretype             = 1                                              # Use exclusion region for channel proteins
[DEBUG  ] |  npbopt               = 0                                              # Use NonLinear PB solver?
[DEBUG  ] |  solvopt              = 1                                              # Select iterative solver
[DEBUG  ] |  accept               = 0.001                                          # Sets the iteration convergence criterion (relative to the initial residue)
[DEBUG  ] |  linit                = 1000                                           # Number of SCF iterations
[DEBUG  ] |  fillratio            = 4.0                                            # Ratio between the longest dimension of the rectangular finite-difference grid and that of the solute
[DEBUG  ] |  scale                = 2.0                                            # 1/scale = grid spacing for the finite difference solver (default = 1/2 Å)
[DEBUG  ] |  nbuffer              = 0.0                                            # Sets how far away (in grid units) the boundary of the finite difference grid is away from the solute surface
[DEBUG  ] |  nfocus               = 2                                              # Electrostatic focusing calculation
[DEBUG  ] |  fscale               = 8                                              # Set the ratio between the coarse and fine grid spacings in an electrostatic focussing calculation
[DEBUG  ] |  npbgrid              = 1                                              # Sets how often the finite-difference grid is regenerated
[DEBUG  ] |  bcopt                = 5                                              # Boundary condition option
[DEBUG  ] |  eneopt               = 2                                              # Compute electrostatic energy and forces
[DEBUG  ] |  frcopt               = 0                                              # Output for computing electrostatic forces
[DEBUG  ] |  scalec               = 0                                              # Option to compute reaction field energy and forces
[DEBUG  ] |  cutfd                = 5.0                                            # Cutoff for finite-difference interactions
[DEBUG  ] |  cutnb                = 0.0                                            # Cutoff for nonbonded interations
[DEBUG  ] |  nsnba                = 1                                              # Sets how often atom-based pairlist is generated
[DEBUG  ] |  decompopt            = 2                                              # Option to select different decomposition schemes when INP = 2
[DEBUG  ] |  use_rmin             = 1                                              # The option to set up van der Waals radii
[DEBUG  ] |  sprob                = 0.557                                          # Solvent probe radius for SASA used to compute the dispersion term
[DEBUG  ] |  vprob                = 1.3                                            # Solvent probe radius for molecular volume (the volume enclosed by SASA)
[DEBUG  ] |  rhow_effect          = 1.129                                          # Effective water density used in the non-polar dispersion term calculation
[DEBUG  ] |  use_sav              = 1                                              # Use molecular volume (the volume enclosed by SASA) for cavity term calculation
[DEBUG  ] |  cavity_surften       = 0.0378                                         # Surface tension
[DEBUG  ] |  cavity_offset        = -0.5692                                        # Offset for nonpolar solvation calc
[DEBUG  ] |  maxsph               = 400                                            # Approximate number of dots to represent the maximum atomic solvent accessible surface
[DEBUG  ] |  maxarcdot            = 1500                                           # Number of dots used to store arc dots per atom
[DEBUG  ] |  npbverb              = 0                                              # Option to turn on verbose mode
[DEBUG  ] |/
[DEBUG  ] |
[DEBUG  ] |# Decomposition namelist variables
[DEBUG  ] |&decomposition
[DEBUG  ] |  idecomp              = 1                                              # Which type of decomposition analysis to do
[DEBUG  ] |  dec_verbose          = 0                                              # Control energy terms are printed to the output
[DEBUG  ] |  print_res            = "within 6"                                     # Which residues to print decomposition data for
[DEBUG  ] |  csv_format           = 1                                              # Write decomposition data in CSV format
[DEBUG  ] |/
[DEBUG  ] |
[DEBUG  ] |--------------------------------------------------------------
[DEBUG  ] 
[INFO   ] Checking mmpbsa.in input file...
[INFO   ] Checking mmpbsa.in input file...Done.

[INFO   ] Checking external programs...
[INFO   ] cpptraj found! Using /home/romi/anaconda3/envs/gmxMMPBSA/bin/cpptraj
[INFO   ] tleap found! Using /home/romi/anaconda3/envs/gmxMMPBSA/bin/tleap
[INFO   ] parmchk2 found! Using /home/romi/anaconda3/envs/gmxMMPBSA/bin/parmchk2
[INFO   ] sander found! Using /home/romi/anaconda3/envs/gmxMMPBSA/bin/sander
[INFO   ] Using GROMACS version > 5.x.x!
[INFO   ] gmx found! Using /usr/local/gromacs/bin/gmx
[INFO   ] Checking external programs...Done.

[INFO   ] Building AMBER topologies from GROMACS files...
[INFO   ] Get PDB files from GROMACS structures files...
[INFO   ] Making gmx_MMPBSA index for complex...
[DEBUG  ] Running command: echo name 18 GMXMMPBSA_REC\n name 19 GMXMMPBSA_LIG\n  18 | 19\n q\n | /usr/local/gromacs/bin/gmx make_ndx -n index.ndx -o _GMXMMPBSA_COM_index.ndx
[DEBUG  ]                      :-) GROMACS - gmx make_ndx, 2021.5 (-:
[DEBUG  ] 
[DEBUG  ]                             GROMACS is written by:
[DEBUG  ]      Andrey Alekseenko              Emile Apol              Rossen Apostolov     
[DEBUG  ]          Paul Bauer           Herman J.C. Berendsen           Par Bjelkmar       
[DEBUG  ]        Christian Blau           Viacheslav Bolnykh             Kevin Boyd        
[DEBUG  ]      Aldert van Buuren           Rudi van Drunen             Anton Feenstra      
[DEBUG  ]     Gilles Gouaillardet             Alan Gray               Gerrit Groenhof      
[DEBUG  ]        Anca Hamuraru            Vincent Hindriksen          M. Eric Irrgang      
[DEBUG  ]       Aleksei Iupinov           Christoph Junghans             Joe Jordan        
[DEBUG  ]     Dimitrios Karkoulis            Peter Kasson                Jiri Kraus        
[DEBUG  ]       Carsten Kutzner              Per Larsson              Justin A. Lemkul     
[DEBUG  ]        Viveca Lindahl            Magnus Lundborg             Erik Marklund       
[DEBUG  ]         Pascal Merz             Pieter Meulenhoff            Teemu Murtola       
[DEBUG  ]         Szilard Pall               Sander Pronk              Roland Schulz       
[DEBUG  ]        Michael Shirts            Alexey Shvetsov             Alfons Sijbers      
[DEBUG  ]        Peter Tieleman              Jon Vincent              Teemu Virolainen     
[DEBUG  ]      Christian Wennberg            Maarten Wolf              Artem Zhmurov       
[DEBUG  ]                            and the project leaders:
[DEBUG  ]         Mark Abraham, Berk Hess, Erik Lindahl, and David van der Spoel
[DEBUG  ] 
[DEBUG  ] Copyright (c) 1991-2000, University of Groningen, The Netherlands.
[DEBUG  ] Copyright (c) 2001-2019, The GROMACS development team at
[DEBUG  ] Uppsala University, Stockholm University and
[DEBUG  ] the Royal Institute of Technology, Sweden.
[DEBUG  ] check out http://www.gromacs.org for more information.
[DEBUG  ] 
[DEBUG  ] GROMACS is free software; you can redistribute it and/or modify it
[DEBUG  ] under the terms of the GNU Lesser General Public License
[DEBUG  ] as published by the Free Software Foundation; either version 2.1
[DEBUG  ] of the License, or (at your option) any later version.
[DEBUG  ] 
[DEBUG  ] GROMACS:      gmx make_ndx, version 2021.5
[DEBUG  ] Executable:   /usr/local/gromacs/bin/gmx
[DEBUG  ] Data prefix:  /usr/local/gromacs
[DEBUG  ] Working dir:  /media/romi/My Book Duo/completed-simulation-oracle/auto-immune/rs2/5mv4-rs2/312k/free-energy
[DEBUG  ] Command line:
[DEBUG  ]   gmx make_ndx -n index.ndx -o _GMXMMPBSA_COM_index.ndx
[DEBUG  ] 
[DEBUG  ] 
[DEBUG  ] GROMACS reminds you: "When It Starts to Start It'll Never Stop" (Magnapop)
[DEBUG  ] 
[DEBUG  ] Going to read 1 old index file(s)
[DEBUG  ] Counted atom numbers up to 163828 in index file
[DEBUG  ] 
[DEBUG  ]   0 System              : 163828 atoms
[DEBUG  ]   1 Protein             :  7012 atoms
[DEBUG  ]   2 Protein-H           :  3571 atoms
[DEBUG  ]   3 C-alpha             :   480 atoms
[DEBUG  ]   4 Backbone            :  1440 atoms
[DEBUG  ]   5 MainChain           :  1917 atoms
[DEBUG  ]   6 MainChain+Cb        :  2344 atoms
[DEBUG  ]   7 MainChain+H         :  2361 atoms
[DEBUG  ]   8 SideChain           :  4651 atoms
[DEBUG  ]   9 SideChain-H         :  1654 atoms
[DEBUG  ]  10 Prot-Masses         :  7012 atoms
[DEBUG  ]  11 non-Protein         : 156816 atoms
[DEBUG  ]  12 Water               : 156813 atoms
[DEBUG  ]  13 SOL                 : 156813 atoms
[DEBUG  ]  14 non-Water           :  7015 atoms
[DEBUG  ]  15 Ion                 :     3 atoms
[DEBUG  ]  16 NA                  :     3 atoms
[DEBUG  ]  17 Water_and_ions      : 156816 atoms
[DEBUG  ]  18 Chain_A             :  6507 atoms
[DEBUG  ]  19 Chain_B             :   505 atoms
[DEBUG  ] 
[DEBUG  ]  nr : group      '!': not  'name' nr name   'splitch' nr    Enter: list groups
[DEBUG  ]  'a': atom       '&': and  'del' nr         'splitres' nr   'l': list residues
[DEBUG  ]  't': atom type  '|': or   'keep' nr        'splitat' nr    'h': help
[DEBUG  ]  'r': residue              'res' nr         'chain' char
[DEBUG  ]  "name": group             'case': case sensitive           'q': save and quit
[DEBUG  ]  'ri': residue index
[DEBUG  ] 
[DEBUG  ] > 
[DEBUG  ] 
[DEBUG  ] > 
[DEBUG  ] 
[DEBUG  ] > 
[DEBUG  ] Copied index group 18 'GMXMMPBSA_REC'
[DEBUG  ] Copied index group 19 'GMXMMPBSA_LIG'
[DEBUG  ] Merged two groups with OR: 6507 505 -> 7012
[DEBUG  ] 
[DEBUG  ]  20 GMXMMPBSA_REC_GMXMMPBSA_LIG:  7012 atoms
[DEBUG  ] 
[DEBUG  ] > 
[INFO   ] Normal Complex: Saving group 18_19 in _GMXMMPBSA_COM_index.ndx file as _GMXMMPBSA_COM.pdb
[DEBUG  ] Running command: echo GMXMMPBSA_REC_GMXMMPBSA_LIG | /usr/local/gromacs/bin/gmx trjconv -f dynamic-nopbc.xtc -s md_0_10.tpr -o _GMXMMPBSA_COM.pdb -n _GMXMMPBSA_COM_index.ndx -dump 0
[DEBUG  ]                      :-) GROMACS - gmx trjconv, 2021.5 (-:
[DEBUG  ] 
[DEBUG  ]                             GROMACS is written by:
[DEBUG  ]      Andrey Alekseenko              Emile Apol              Rossen Apostolov     
[DEBUG  ]          Paul Bauer           Herman J.C. Berendsen           Par Bjelkmar       
[DEBUG  ]        Christian Blau           Viacheslav Bolnykh             Kevin Boyd        
[DEBUG  ]      Aldert van Buuren           Rudi van Drunen             Anton Feenstra      
[DEBUG  ]     Gilles Gouaillardet             Alan Gray               Gerrit Groenhof      
[DEBUG  ]        Anca Hamuraru            Vincent Hindriksen          M. Eric Irrgang      
[DEBUG  ]       Aleksei Iupinov           Christoph Junghans             Joe Jordan        
[DEBUG  ]     Dimitrios Karkoulis            Peter Kasson                Jiri Kraus        
[DEBUG  ]       Carsten Kutzner              Per Larsson              Justin A. Lemkul     
[DEBUG  ]        Viveca Lindahl            Magnus Lundborg             Erik Marklund       
[DEBUG  ]         Pascal Merz             Pieter Meulenhoff            Teemu Murtola       
[DEBUG  ]         Szilard Pall               Sander Pronk              Roland Schulz       
[DEBUG  ]        Michael Shirts            Alexey Shvetsov             Alfons Sijbers      
[DEBUG  ]        Peter Tieleman              Jon Vincent              Teemu Virolainen     
[DEBUG  ]      Christian Wennberg            Maarten Wolf              Artem Zhmurov       
[DEBUG  ]                            and the project leaders:
[DEBUG  ]         Mark Abraham, Berk Hess, Erik Lindahl, and David van der Spoel
[DEBUG  ] 
[DEBUG  ] Copyright (c) 1991-2000, University of Groningen, The Netherlands.
[DEBUG  ] Copyright (c) 2001-2019, The GROMACS development team at
[DEBUG  ] Uppsala University, Stockholm University and
[DEBUG  ] the Royal Institute of Technology, Sweden.
[DEBUG  ] check out http://www.gromacs.org for more information.
[DEBUG  ] 
[DEBUG  ] GROMACS is free software; you can redistribute it and/or modify it
[DEBUG  ] under the terms of the GNU Lesser General Public License
[DEBUG  ] as published by the Free Software Foundation; either version 2.1
[DEBUG  ] of the License, or (at your option) any later version.
[DEBUG  ] 
[DEBUG  ] GROMACS:      gmx trjconv, version 2021.5
[DEBUG  ] Executable:   /usr/local/gromacs/bin/gmx
[DEBUG  ] Data prefix:  /usr/local/gromacs
[DEBUG  ] Working dir:  /media/romi/My Book Duo/completed-simulation-oracle/auto-immune/rs2/5mv4-rs2/312k/free-energy
[DEBUG  ] Command line:
[DEBUG  ]   gmx trjconv -f dynamic-nopbc.xtc -s md_0_10.tpr -o _GMXMMPBSA_COM.pdb -n _GMXMMPBSA_COM_index.ndx -dump 0
[DEBUG  ] 
[DEBUG  ] Will write pdb: Protein data bank file
[DEBUG  ] Reading file md_0_10.tpr, VERSION 2019.6 (single precision)
[DEBUG  ] Reading file md_0_10.tpr, VERSION 2019.6 (single precision)
[DEBUG  ] Group     0 (         System) has 163828 elements
[DEBUG  ] Group     1 (        Protein) has  7012 elements
[DEBUG  ] Group     2 (      Protein-H) has  3571 elements
[DEBUG  ] Group     3 (        C-alpha) has   480 elements
[DEBUG  ] Group     4 (       Backbone) has  1440 elements
[DEBUG  ] Group     5 (      MainChain) has  1917 elements
[DEBUG  ] Group     6 (   MainChain+Cb) has  2344 elements
[DEBUG  ] Group     7 (    MainChain+H) has  2361 elements
[DEBUG  ] Group     8 (      SideChain) has  4651 elements
[DEBUG  ] Group     9 (    SideChain-H) has  1654 elements
[DEBUG  ] Group    10 (    Prot-Masses) has  7012 elements
[DEBUG  ] Group    11 (    non-Protein) has 156816 elements
[DEBUG  ] Group    12 (          Water) has 156813 elements
[DEBUG  ] Group    13 (            SOL) has 156813 elements
[DEBUG  ] Group    14 (      non-Water) has  7015 elements
[DEBUG  ] Group    15 (            Ion) has     3 elements
[DEBUG  ] Group    16 (             NA) has     3 elements
[DEBUG  ] Group    17 ( Water_and_ions) has 156816 elements
[DEBUG  ] Group    18 (  GMXMMPBSA_REC) has  6507 elements
[DEBUG  ] Group    19 (  GMXMMPBSA_LIG) has   505 elements
[DEBUG  ] Group    20 (GMXMMPBSA_REC_GMXMMPBSA_LIG) has  7012 elements
[DEBUG  ] Select a group: Reading frame       0 time    0.000   
[DEBUG  ] Precision of dynamic-nopbc.xtc is 0.001 (nm)
[DEBUG  ] Reading frame       1 time   10.000   Reading frame       0 time    0.000   
[DEBUG  ] Dumping frame at t= 0 ps
[DEBUG  ]  ->  frame      0 time    0.000      Reading frame       1 time   10.000   
[DEBUG  ] 
[DEBUG  ] GROMACS reminds you: "When It Starts to Start It'll Never Stop" (Magnapop)
[DEBUG  ] 
[DEBUG  ] Note that major changes are planned in future for trjconv, to improve usability and utility.
[DEBUG  ] Select group for output
[DEBUG  ] Selected 20: 'GMXMMPBSA_REC_GMXMMPBSA_LIG'
[INFO   ] No receptor structure file was defined. Using ST approach...
[INFO   ] Using receptor structure from complex to generate AMBER topology
[INFO   ] Normal Receptor: Saving group 18 in _GMXMMPBSA_COM_index.ndx file as _GMXMMPBSA_REC.pdb
[DEBUG  ] Running command: echo 18 | /usr/local/gromacs/bin/gmx trjconv -f dynamic-nopbc.xtc -s md_0_10.tpr -o _GMXMMPBSA_REC.pdb -n _GMXMMPBSA_COM_index.ndx -dump 0
[DEBUG  ]                      :-) GROMACS - gmx trjconv, 2021.5 (-:
[DEBUG  ] 
[DEBUG  ]                             GROMACS is written by:
[DEBUG  ]      Andrey Alekseenko              Emile Apol              Rossen Apostolov     
[DEBUG  ]          Paul Bauer           Herman J.C. Berendsen           Par Bjelkmar       
[DEBUG  ]        Christian Blau           Viacheslav Bolnykh             Kevin Boyd        
[DEBUG  ]      Aldert van Buuren           Rudi van Drunen             Anton Feenstra      
[DEBUG  ]     Gilles Gouaillardet             Alan Gray               Gerrit Groenhof      
[DEBUG  ]        Anca Hamuraru            Vincent Hindriksen          M. Eric Irrgang      
[DEBUG  ]       Aleksei Iupinov           Christoph Junghans             Joe Jordan        
[DEBUG  ]     Dimitrios Karkoulis            Peter Kasson                Jiri Kraus        
[DEBUG  ]       Carsten Kutzner              Per Larsson              Justin A. Lemkul     
[DEBUG  ]        Viveca Lindahl            Magnus Lundborg             Erik Marklund       
[DEBUG  ]         Pascal Merz             Pieter Meulenhoff            Teemu Murtola       
[DEBUG  ]         Szilard Pall               Sander Pronk              Roland Schulz       
[DEBUG  ]        Michael Shirts            Alexey Shvetsov             Alfons Sijbers      
[DEBUG  ]        Peter Tieleman              Jon Vincent              Teemu Virolainen     
[DEBUG  ]      Christian Wennberg            Maarten Wolf              Artem Zhmurov       
[DEBUG  ]                            and the project leaders:
[DEBUG  ]         Mark Abraham, Berk Hess, Erik Lindahl, and David van der Spoel
[DEBUG  ] 
[DEBUG  ] Copyright (c) 1991-2000, University of Groningen, The Netherlands.
[DEBUG  ] Copyright (c) 2001-2019, The GROMACS development team at
[DEBUG  ] Uppsala University, Stockholm University and
[DEBUG  ] the Royal Institute of Technology, Sweden.
[DEBUG  ] check out http://www.gromacs.org for more information.
[DEBUG  ] 
[DEBUG  ] GROMACS is free software; you can redistribute it and/or modify it
[DEBUG  ] under the terms of the GNU Lesser General Public License
[DEBUG  ] as published by the Free Software Foundation; either version 2.1
[DEBUG  ] of the License, or (at your option) any later version.
[DEBUG  ] 
[DEBUG  ] GROMACS:      gmx trjconv, version 2021.5
[DEBUG  ] Executable:   /usr/local/gromacs/bin/gmx
[DEBUG  ] Data prefix:  /usr/local/gromacs
[DEBUG  ] Working dir:  /media/romi/My Book Duo/completed-simulation-oracle/auto-immune/rs2/5mv4-rs2/312k/free-energy
[DEBUG  ] Command line:
[DEBUG  ]   gmx trjconv -f dynamic-nopbc.xtc -s md_0_10.tpr -o _GMXMMPBSA_REC.pdb -n _GMXMMPBSA_COM_index.ndx -dump 0
[DEBUG  ] 
[DEBUG  ] Will write pdb: Protein data bank file
[DEBUG  ] Reading file md_0_10.tpr, VERSION 2019.6 (single precision)
[DEBUG  ] Reading file md_0_10.tpr, VERSION 2019.6 (single precision)
[DEBUG  ] Group     0 (         System) has 163828 elements
[DEBUG  ] Group     1 (        Protein) has  7012 elements
[DEBUG  ] Group     2 (      Protein-H) has  3571 elements
[DEBUG  ] Group     3 (        C-alpha) has   480 elements
[DEBUG  ] Group     4 (       Backbone) has  1440 elements
[DEBUG  ] Group     5 (      MainChain) has  1917 elements
[DEBUG  ] Group     6 (   MainChain+Cb) has  2344 elements
[DEBUG  ] Group     7 (    MainChain+H) has  2361 elements
[DEBUG  ] Group     8 (      SideChain) has  4651 elements
[DEBUG  ] Group     9 (    SideChain-H) has  1654 elements
[DEBUG  ] Group    10 (    Prot-Masses) has  7012 elements
[DEBUG  ] Group    11 (    non-Protein) has 156816 elements
[DEBUG  ] Group    12 (          Water) has 156813 elements
[DEBUG  ] Group    13 (            SOL) has 156813 elements
[DEBUG  ] Group    14 (      non-Water) has  7015 elements
[DEBUG  ] Group    15 (            Ion) has     3 elements
[DEBUG  ] Group    16 (             NA) has     3 elements
[DEBUG  ] Group    17 ( Water_and_ions) has 156816 elements
[DEBUG  ] Group    18 (  GMXMMPBSA_REC) has  6507 elements
[DEBUG  ] Group    19 (  GMXMMPBSA_LIG) has   505 elements
[DEBUG  ] Group    20 (GMXMMPBSA_REC_GMXMMPBSA_LIG) has  7012 elements
[DEBUG  ] Select a group: Reading frame       0 time    0.000   
[DEBUG  ] Precision of dynamic-nopbc.xtc is 0.001 (nm)
[DEBUG  ] Reading frame       1 time   10.000   Reading frame       0 time    0.000   
[DEBUG  ] Dumping frame at t= 0 ps
[DEBUG  ]  ->  frame      0 time    0.000      Reading frame       1 time   10.000   
[DEBUG  ] 
[DEBUG  ] GROMACS reminds you: "When It Starts to Start It'll Never Stop" (Magnapop)
[DEBUG  ] 
[DEBUG  ] Note that major changes are planned in future for trjconv, to improve usability and utility.
[DEBUG  ] Select group for output
[DEBUG  ] Selected 18: 'GMXMMPBSA_REC'
[INFO   ] No ligand structure file was defined. Using ST approach...
[INFO   ] Using ligand structure from complex to generate AMBER topology
[INFO   ] Normal Ligand: Saving group 19 in _GMXMMPBSA_COM_index.ndx file as _GMXMMPBSA_LIG.pdb
[DEBUG  ] Running command: echo 19 | /usr/local/gromacs/bin/gmx trjconv -f dynamic-nopbc.xtc -s md_0_10.tpr -o _GMXMMPBSA_LIG.pdb -n _GMXMMPBSA_COM_index.ndx -dump 0
[DEBUG  ]                      :-) GROMACS - gmx trjconv, 2021.5 (-:
[DEBUG  ] 
[DEBUG  ]                             GROMACS is written by:
[DEBUG  ]      Andrey Alekseenko              Emile Apol              Rossen Apostolov     
[DEBUG  ]          Paul Bauer           Herman J.C. Berendsen           Par Bjelkmar       
[DEBUG  ]        Christian Blau           Viacheslav Bolnykh             Kevin Boyd        
[DEBUG  ]      Aldert van Buuren           Rudi van Drunen             Anton Feenstra      
[DEBUG  ]     Gilles Gouaillardet             Alan Gray               Gerrit Groenhof      
[DEBUG  ]        Anca Hamuraru            Vincent Hindriksen          M. Eric Irrgang      
[DEBUG  ]       Aleksei Iupinov           Christoph Junghans             Joe Jordan        
[DEBUG  ]     Dimitrios Karkoulis            Peter Kasson                Jiri Kraus        
[DEBUG  ]       Carsten Kutzner              Per Larsson              Justin A. Lemkul     
[DEBUG  ]        Viveca Lindahl            Magnus Lundborg             Erik Marklund       
[DEBUG  ]         Pascal Merz             Pieter Meulenhoff            Teemu Murtola       
[DEBUG  ]         Szilard Pall               Sander Pronk              Roland Schulz       
[DEBUG  ]        Michael Shirts            Alexey Shvetsov             Alfons Sijbers      
[DEBUG  ]        Peter Tieleman              Jon Vincent              Teemu Virolainen     
[DEBUG  ]      Christian Wennberg            Maarten Wolf              Artem Zhmurov       
[DEBUG  ]                            and the project leaders:
[DEBUG  ]         Mark Abraham, Berk Hess, Erik Lindahl, and David van der Spoel
[DEBUG  ] 
[DEBUG  ] Copyright (c) 1991-2000, University of Groningen, The Netherlands.
[DEBUG  ] Copyright (c) 2001-2019, The GROMACS development team at
[DEBUG  ] Uppsala University, Stockholm University and
[DEBUG  ] the Royal Institute of Technology, Sweden.
[DEBUG  ] check out http://www.gromacs.org for more information.
[DEBUG  ] 
[DEBUG  ] GROMACS is free software; you can redistribute it and/or modify it
[DEBUG  ] under the terms of the GNU Lesser General Public License
[DEBUG  ] as published by the Free Software Foundation; either version 2.1
[DEBUG  ] of the License, or (at your option) any later version.
[DEBUG  ] 
[DEBUG  ] GROMACS:      gmx trjconv, version 2021.5
[DEBUG  ] Executable:   /usr/local/gromacs/bin/gmx
[DEBUG  ] Data prefix:  /usr/local/gromacs
[DEBUG  ] Working dir:  /media/romi/My Book Duo/completed-simulation-oracle/auto-immune/rs2/5mv4-rs2/312k/free-energy
[DEBUG  ] Command line:
[DEBUG  ]   gmx trjconv -f dynamic-nopbc.xtc -s md_0_10.tpr -o _GMXMMPBSA_LIG.pdb -n _GMXMMPBSA_COM_index.ndx -dump 0
[DEBUG  ] 
[DEBUG  ] Will write pdb: Protein data bank file
[DEBUG  ] Reading file md_0_10.tpr, VERSION 2019.6 (single precision)
[DEBUG  ] Reading file md_0_10.tpr, VERSION 2019.6 (single precision)
[DEBUG  ] Group     0 (         System) has 163828 elements
[DEBUG  ] Group     1 (        Protein) has  7012 elements
[DEBUG  ] Group     2 (      Protein-H) has  3571 elements
[DEBUG  ] Group     3 (        C-alpha) has   480 elements
[DEBUG  ] Group     4 (       Backbone) has  1440 elements
[DEBUG  ] Group     5 (      MainChain) has  1917 elements
[DEBUG  ] Group     6 (   MainChain+Cb) has  2344 elements
[DEBUG  ] Group     7 (    MainChain+H) has  2361 elements
[DEBUG  ] Group     8 (      SideChain) has  4651 elements
[DEBUG  ] Group     9 (    SideChain-H) has  1654 elements
[DEBUG  ] Group    10 (    Prot-Masses) has  7012 elements
[DEBUG  ] Group    11 (    non-Protein) has 156816 elements
[DEBUG  ] Group    12 (          Water) has 156813 elements
[DEBUG  ] Group    13 (            SOL) has 156813 elements
[DEBUG  ] Group    14 (      non-Water) has  7015 elements
[DEBUG  ] Group    15 (            Ion) has     3 elements
[DEBUG  ] Group    16 (             NA) has     3 elements
[DEBUG  ] Group    17 ( Water_and_ions) has 156816 elements
[DEBUG  ] Group    18 (  GMXMMPBSA_REC) has  6507 elements
[DEBUG  ] Group    19 (  GMXMMPBSA_LIG) has   505 elements
[DEBUG  ] Group    20 (GMXMMPBSA_REC_GMXMMPBSA_LIG) has  7012 elements
[DEBUG  ] Select a group: Reading frame       0 time    0.000   
[DEBUG  ] Precision of dynamic-nopbc.xtc is 0.001 (nm)
[DEBUG  ] Reading frame       1 time   10.000   Reading frame       0 time    0.000   
[DEBUG  ] Dumping frame at t= 0 ps
[DEBUG  ]  ->  frame      0 time    0.000      Reading frame       1 time   10.000   
[DEBUG  ] 
[DEBUG  ] GROMACS reminds you: "When It Starts to Start It'll Never Stop" (Magnapop)
[DEBUG  ] 
[DEBUG  ] Note that major changes are planned in future for trjconv, to improve usability and utility.
[DEBUG  ] Select group for output
[DEBUG  ] Selected 19: 'GMXMMPBSA_LIG'
[INFO   ] Checking the structures consistency...
[INFO   ] 
[INFO   ] Using topology conversion. Setting radiopt = 0...
[INFO   ] Building Normal Complex Amber topology...
[INFO   ] Detected CHARMM force field topology format...
[INFO   ] Assigning PBRadii mbondi2 to Complex...
[INFO   ] Writing Normal Complex AMBER topology...
[INFO   ] No Receptor topology file was defined. Using ST approach...
[INFO   ] Building AMBER Receptor topology from Complex...
[INFO   ] Assigning PBRadii mbondi2 to Receptor...
[INFO   ] Writing Normal Receptor AMBER topology...
[INFO   ] No Ligand topology file was defined. Using ST approach...
[INFO   ] Building AMBER Ligand topology from Complex...
[INFO   ] Assigning PBRadii mbondi2 to Ligand...
[INFO   ] Writing Normal Ligand AMBER topology...
[INFO   ] Selecting residues by distance (6 Å) between receptor and ligand for decomposition analysis...
[INFO   ] Selected 59 residues:
R:A:PHE:27  R:A:THR:28  R:A:SER:30  R:A:ASP:31  R:A:ALA:32  R:A:TRP:33  R:A:GLU:50  R:A:ARG:52  R:A:ASN:53  R:A:VAL:55
R:A:ASN:56  R:A:GLY:100 R:A:LEU:101 R:A:THR:102 R:A:PHE:103 R:A:ASP:104 R:A:TYR:105 R:A:TRP:106 R:B:ASN:249 R:B:TYR:250
R:B:GLY:251 R:B:ILE:252 R:B:SER:254 R:B:ASN:256 R:B:TRP:257 R:B:PHE:258 R:B:PHE:268 R:B:LEU:269 R:B:ILE:270 R:B:TYR:271
R:B:ALA:272 R:B:LYS:275 R:B:GLN:276 R:B:GLY:277 R:B:SER:278 R:B:GLN:311 R:B:SER:313 R:B:LYS:314 R:B:GLY:315 R:B:VAL:316
R:B:TYR:318 L:C:ALA:456 L:C:GLY:457 L:C:GLU:458 L:C:ARG:459 L:C:GLY:460 L:C:GLU:461 L:C:THR:462 L:C:GLY:463 L:C:PRO:464
L:C:GLY:465 L:C:PRO:466 L:C:ALA:467 L:C:GLY:468 L:C:PHE:469 L:C:ALA:470 L:C:GLY:471 L:C:PRO:472 L:C:PRO:473

[INFO   ] Cleaning normal complex trajectories...
[DEBUG  ] Running command: echo GMXMMPBSA_REC_GMXMMPBSA_LIG | /usr/local/gromacs/bin/gmx trjconv -f dynamic-nopbc.xtc -s md_0_10.tpr -o COM_traj_0.xtc -n _GMXMMPBSA_COM_index.ndx
[DEBUG  ]                      :-) GROMACS - gmx trjconv, 2021.5 (-:
[DEBUG  ] 
[DEBUG  ]                             GROMACS is written by:
[DEBUG  ]      Andrey Alekseenko              Emile Apol              Rossen Apostolov     
[DEBUG  ]          Paul Bauer           Herman J.C. Berendsen           Par Bjelkmar       
[DEBUG  ]        Christian Blau           Viacheslav Bolnykh             Kevin Boyd        
[DEBUG  ]      Aldert van Buuren           Rudi van Drunen             Anton Feenstra      
[DEBUG  ]     Gilles Gouaillardet             Alan Gray               Gerrit Groenhof      
[DEBUG  ]        Anca Hamuraru            Vincent Hindriksen          M. Eric Irrgang      
[DEBUG  ]       Aleksei Iupinov           Christoph Junghans             Joe Jordan        
[DEBUG  ]     Dimitrios Karkoulis            Peter Kasson                Jiri Kraus        
[DEBUG  ]       Carsten Kutzner              Per Larsson              Justin A. Lemkul     
[DEBUG  ]        Viveca Lindahl            Magnus Lundborg             Erik Marklund       
[DEBUG  ]         Pascal Merz             Pieter Meulenhoff            Teemu Murtola       
[DEBUG  ]         Szilard Pall               Sander Pronk              Roland Schulz       
[DEBUG  ]        Michael Shirts            Alexey Shvetsov             Alfons Sijbers      
[DEBUG  ]        Peter Tieleman              Jon Vincent              Teemu Virolainen     
[DEBUG  ]      Christian Wennberg            Maarten Wolf              Artem Zhmurov       
[DEBUG  ]                            and the project leaders:
[DEBUG  ]         Mark Abraham, Berk Hess, Erik Lindahl, and David van der Spoel
[DEBUG  ] 
[DEBUG  ] Copyright (c) 1991-2000, University of Groningen, The Netherlands.
[DEBUG  ] Copyright (c) 2001-2019, The GROMACS development team at
[DEBUG  ] Uppsala University, Stockholm University and
[DEBUG  ] the Royal Institute of Technology, Sweden.
[DEBUG  ] check out http://www.gromacs.org for more information.
[DEBUG  ] 
[DEBUG  ] GROMACS is free software; you can redistribute it and/or modify it
[DEBUG  ] under the terms of the GNU Lesser General Public License
[DEBUG  ] as published by the Free Software Foundation; either version 2.1
[DEBUG  ] of the License, or (at your option) any later version.
[DEBUG  ] 
[DEBUG  ] GROMACS:      gmx trjconv, version 2021.5
[DEBUG  ] Executable:   /usr/local/gromacs/bin/gmx
[DEBUG  ] Data prefix:  /usr/local/gromacs
[DEBUG  ] Working dir:  /media/romi/My Book Duo/completed-simulation-oracle/auto-immune/rs2/5mv4-rs2/312k/free-energy
[DEBUG  ] Command line:
[DEBUG  ]   gmx trjconv -f dynamic-nopbc.xtc -s md_0_10.tpr -o COM_traj_0.xtc -n _GMXMMPBSA_COM_index.ndx
[DEBUG  ] 
[DEBUG  ] Will write xtc: Compressed trajectory (portable xdr format): xtc
[DEBUG  ] Reading file md_0_10.tpr, VERSION 2019.6 (single precision)
[DEBUG  ] Reading file md_0_10.tpr, VERSION 2019.6 (single precision)
[DEBUG  ] Group     0 (         System) has 163828 elements
[DEBUG  ] Group     1 (        Protein) has  7012 elements
[DEBUG  ] Group     2 (      Protein-H) has  3571 elements
[DEBUG  ] Group     3 (        C-alpha) has   480 elements
[DEBUG  ] Group     4 (       Backbone) has  1440 elements
[DEBUG  ] Group     5 (      MainChain) has  1917 elements
[DEBUG  ] Group     6 (   MainChain+Cb) has  2344 elements
[DEBUG  ] Group     7 (    MainChain+H) has  2361 elements
[DEBUG  ] Group     8 (      SideChain) has  4651 elements
[DEBUG  ] Group     9 (    SideChain-H) has  1654 elements
[DEBUG  ] Group    10 (    Prot-Masses) has  7012 elements
[DEBUG  ] Group    11 (    non-Protein) has 156816 elements
[DEBUG  ] Group    12 (          Water) has 156813 elements
[DEBUG  ] Group    13 (            SOL) has 156813 elements
[DEBUG  ] Group    14 (      non-Water) has  7015 elements
[DEBUG  ] Group    15 (            Ion) has     3 elements
[DEBUG  ] Group    16 (             NA) has     3 elements
[DEBUG  ] Group    17 ( Water_and_ions) has 156816 elements
[DEBUG  ] Group    18 (  GMXMMPBSA_REC) has  6507 elements
[DEBUG  ] Group    19 (  GMXMMPBSA_LIG) has   505 elements
[DEBUG  ] Group    20 (GMXMMPBSA_REC_GMXMMPBSA_LIG) has  7012 elements
[DEBUG  ] Select a group: Reading frame       0 time    0.000   
[DEBUG  ] Precision of dynamic-nopbc.xtc is 0.001 (nm)
[DEBUG  ] Using output precision of 0.001 (nm)
[DEBUG  ]  ->  frame      0 time    0.000      Reading frame       1 time   10.000    ->  frame      1 time   10.000      Reading frame       2 time   20.000    ->  frame      2 time   20.000      Reading frame       3 time   30.000    ->  frame      3 time   30.000      Reading frame       4 time   40.000    ->  frame      4 time   40.000      Reading frame       5 time   50.000    ->  frame      5 time   50.000      Reading frame       6 time   60.000    ->  frame      6 time   60.000      Reading frame       7 time   70.000    ->  frame      7 time   70.000      Reading frame       8 time   80.000    ->  frame      8 time   80.000      Reading frame       9 time   90.000    ->  frame      9 time   90.000      Reading frame      10 time  100.000    ->  frame     10 time  100.000      Reading frame      11 time  110.000   Reading frame      12 time  120.000   Reading frame      13 time  130.000   Reading frame      14 time  140.000   Reading frame      15 time  150.000   Reading frame      16 time  160.000   Reading frame      17 time  170.000   Reading frame      18 time  180.000   Reading frame      19 time  190.000   Reading frame      20 time  200.000    ->  frame     20 time  200.000      Reading frame      30 time  300.000    ->  frame     30 time  300.000      Reading frame      40 time  400.000    ->  frame     40 time  400.000      Reading frame      50 time  500.000    ->  frame     50 time  500.000      Reading frame      60 time  600.000    ->  frame     60 time  600.000      Reading frame      70 time  700.000    ->  frame     70 time  700.000      Reading frame      80 time  800.000    ->  frame     80 time  800.000      Reading frame      90 time  900.000    ->  frame     90 time  900.000      Reading frame     100 time 1000.000    ->  frame    100 time 1000.000      Reading frame     110 time 1100.000    ->  frame    110 time 1100.000      Reading frame     120 time 1200.000    ->  frame    120 time 1200.000      Reading frame     130 time 1300.000    ->  frame    130 time 1300.000      Reading frame     140 time 1400.000    ->  frame    140 time 1400.000      Reading frame     150 time 1500.000    ->  frame    150 time 1500.000      Reading frame     160 time 1600.000    ->  frame    160 time 1600.000      Reading frame     170 time 1700.000    ->  frame    170 time 1700.000      Reading frame     180 time 1800.000    ->  frame    180 time 1800.000      Reading frame     190 time 1900.000    ->  frame    190 time 1900.000      Reading frame     200 time 2000.000    ->  frame    200 time 2000.000       ->  frame    210 time 2100.000       ->  frame    220 time 2200.000       ->  frame    230 time 2300.000       ->  frame    240 time 2400.000       ->  frame    250 time 2500.000       ->  frame    260 time 2600.000       ->  frame    270 time 2700.000       ->  frame    280 time 2800.000       ->  frame    290 time 2900.000      Reading frame     300 time 3000.000    ->  frame    300 time 3000.000       ->  frame    310 time 3100.000       ->  frame    320 time 3200.000       ->  frame    330 time 3300.000       ->  frame    340 time 3400.000       ->  frame    350 time 3500.000       ->  frame    360 time 3600.000       ->  frame    370 time 3700.000       ->  frame    380 time 3800.000       ->  frame    390 time 3900.000      Reading frame     400 time 4000.000    ->  frame    400 time 4000.000       ->  frame    410 time 4100.000       ->  frame    420 time 4200.000       ->  frame    430 time 4300.000       ->  frame    440 time 4400.000       ->  frame    450 time 4500.000       ->  frame    460 time 4600.000       ->  frame    470 time 4700.000       ->  frame    480 time 4800.000       ->  frame    490 time 4900.000      Reading frame     500 time 5000.000    ->  frame    500 time 5000.000       ->  frame    510 time 5100.000       ->  frame    520 time 5200.000       ->  frame    530 time 5300.000       ->  frame    540 time 5400.000       ->  frame    550 time 5500.000       ->  frame    560 time 5600.000       ->  frame    570 time 5700.000       ->  frame    580 time 5800.000       ->  frame    590 time 5900.000      Reading frame     600 time 6000.000    ->  frame    600 time 6000.000       ->  frame    610 time 6100.000       ->  frame    620 time 6200.000       ->  frame    630 time 6300.000       ->  frame    640 time 6400.000       ->  frame    650 time 6500.000       ->  frame    660 time 6600.000       ->  frame    670 time 6700.000       ->  frame    680 time 6800.000       ->  frame    690 time 6900.000      Reading frame     700 time 7000.000    ->  frame    700 time 7000.000       ->  frame    710 time 7100.000       ->  frame    720 time 7200.000       ->  frame    730 time 7300.000       ->  frame    740 time 7400.000       ->  frame    750 time 7500.000       ->  frame    760 time 7600.000       ->  frame    770 time 7700.000       ->  frame    780 time 7800.000       ->  frame    790 time 7900.000      Reading frame     800 time 8000.000    ->  frame    800 time 8000.000       ->  frame    810 time 8100.000       ->  frame    820 time 8200.000       ->  frame    830 time 8300.000       ->  frame    840 time 8400.000       ->  frame    850 time 8500.000       ->  frame    860 time 8600.000       ->  frame    870 time 8700.000       ->  frame    880 time 8800.000       ->  frame    890 time 8900.000      Reading frame     900 time 9000.000    ->  frame    900 time 9000.000       ->  frame    910 time 9100.000       ->  frame    920 time 9200.000       ->  frame    930 time 9300.000       ->  frame    940 time 9400.000       ->  frame    950 time 9500.000       ->  frame    960 time 9600.000       ->  frame    970 time 9700.000       ->  frame    980 time 9800.000       ->  frame    990 time 9900.000      Reading frame    1000 time 10000.000    ->  frame   1000 time 10000.000       ->  frame   1010 time 10100.000       ->  frame   1020 time 10200.000       ->  frame   1030 time 10300.000       ->  frame   1040 time 10400.000       ->  frame   1050 time 10500.000       ->  frame   1060 time 10600.000       ->  frame   1070 time 10700.000       ->  frame   1080 time 10800.000       ->  frame   1090 time 10900.000      Reading frame    1100 time 11000.000    ->  frame   1100 time 11000.000       ->  frame   1110 time 11100.000       ->  frame   1120 time 11200.000       ->  frame   1130 time 11300.000       ->  frame   1140 time 11400.000       ->  frame   1150 time 11500.000       ->  frame   1160 time 11600.000       ->  frame   1170 time 11700.000       ->  frame   1180 time 11800.000       ->  frame   1190 time 11900.000      Reading frame    1200 time 12000.000    ->  frame   1200 time 12000.000       ->  frame   1210 time 12100.000       ->  frame   1220 time 12200.000       ->  frame   1230 time 12300.000       ->  frame   1240 time 12400.000       ->  frame   1250 time 12500.000       ->  frame   1260 time 12600.000       ->  frame   1270 time 12700.000       ->  frame   1280 time 12800.000       ->  frame   1290 time 12900.000      Reading frame    1300 time 13000.000    ->  frame   1300 time 13000.000       ->  frame   1310 time 13100.000       ->  frame   1320 time 13200.000       ->  frame   1330 time 13300.000       ->  frame   1340 time 13400.000       ->  frame   1350 time 13500.000       ->  frame   1360 time 13600.000       ->  frame   1370 time 13700.000       ->  frame   1380 time 13800.000       ->  frame   1390 time 13900.000      Reading frame    1400 time 14000.000    ->  frame   1400 time 14000.000       ->  frame   1410 time 14100.000       ->  frame   1420 time 14200.000       ->  frame   1430 time 14300.000       ->  frame   1440 time 14400.000       ->  frame   1450 time 14500.000       ->  frame   1460 time 14600.000       ->  frame   1470 time 14700.000       ->  frame   1480 time 14800.000       ->  frame   1490 time 14900.000      Reading frame    1500 time 15000.000    ->  frame   1500 time 15000.000       ->  frame   1510 time 15100.000       ->  frame   1520 time 15200.000       ->  frame   1530 time 15300.000       ->  frame   1540 time 15400.000       ->  frame   1550 time 15500.000       ->  frame   1560 time 15600.000       ->  frame   1570 time 15700.000       ->  frame   1580 time 15800.000       ->  frame   1590 time 15900.000      Reading frame    1600 time 16000.000    ->  frame   1600 time 16000.000       ->  frame   1610 time 16100.000       ->  frame   1620 time 16200.000       ->  frame   1630 time 16300.000       ->  frame   1640 time 16400.000       ->  frame   1650 time 16500.000       ->  frame   1660 time 16600.000       ->  frame   1670 time 16700.000       ->  frame   1680 time 16800.000       ->  frame   1690 time 16900.000      Reading frame    1700 time 17000.000    ->  frame   1700 time 17000.000       ->  frame   1710 time 17100.000       ->  frame   1720 time 17200.000       ->  frame   1730 time 17300.000       ->  frame   1740 time 17400.000       ->  frame   1750 time 17500.000       ->  frame   1760 time 17600.000       ->  frame   1770 time 17700.000       ->  frame   1780 time 17800.000       ->  frame   1790 time 17900.000      Reading frame    1800 time 18000.000    ->  frame   1800 time 18000.000       ->  frame   1810 time 18100.000       ->  frame   1820 time 18200.000       ->  frame   1830 time 18300.000       ->  frame   1840 time 18400.000       ->  frame   1850 time 18500.000       ->  frame   1860 time 18600.000       ->  frame   1870 time 18700.000       ->  frame   1880 time 18800.000       ->  frame   1890 time 18900.000      Reading frame    1900 time 19000.000    ->  frame   1900 time 19000.000       ->  frame   1910 time 19100.000       ->  frame   1920 time 19200.000       ->  frame   1930 time 19300.000       ->  frame   1940 time 19400.000       ->  frame   1950 time 19500.000       ->  frame   1960 time 19600.000       ->  frame   1970 time 19700.000       ->  frame   1980 time 19800.000       ->  frame   1990 time 19900.000      Reading frame    2000 time 20000.000    ->  frame   2000 time 20000.000       ->  frame   2010 time 20100.000       ->  frame   2020 time 20200.000       ->  frame   2030 time 20300.000       ->  frame   2040 time 20400.000       ->  frame   2050 time 20500.000       ->  frame   2060 time 20600.000       ->  frame   2070 time 20700.000       ->  frame   2080 time 20800.000       ->  frame   2090 time 20900.000       ->  frame   2100 time 21000.000       ->  frame   2110 time 21100.000       ->  frame   2120 time 21200.000       ->  frame   2130 time 21300.000       ->  frame   2140 time 21400.000       ->  frame   2150 time 21500.000       ->  frame   2160 time 21600.000       ->  frame   2170 time 21700.000       ->  frame   2180 time 21800.000       ->  frame   2190 time 21900.000       ->  frame   2200 time 22000.000       ->  frame   2210 time 22100.000       ->  frame   2220 time 22200.000       ->  frame   2230 time 22300.000       ->  frame   2240 time 22400.000       ->  frame   2250 time 22500.000       ->  frame   2260 time 22600.000       ->  frame   2270 time 22700.000       ->  frame   2280 time 22800.000       ->  frame   2290 time 22900.000       ->  frame   2300 time 23000.000       ->  frame   2310 time 23100.000       ->  frame   2320 time 23200.000       ->  frame   2330 time 23300.000       ->  frame   2340 time 23400.000       ->  frame   2350 time 23500.000       ->  frame   2360 time 23600.000       ->  frame   2370 time 23700.000       ->  frame   2380 time 23800.000       ->  frame   2390 time 23900.000       ->  frame   2400 time 24000.000       ->  frame   2410 time 24100.000       ->  frame   2420 time 24200.000       ->  frame   2430 time 24300.000       ->  frame   2440 time 24400.000       ->  frame   2450 time 24500.000       ->  frame   2460 time 24600.000       ->  frame   2470 time 24700.000       ->  frame   2480 time 24800.000       ->  frame   2490 time 24900.000       ->  frame   2500 time 25000.000       ->  frame   2510 time 25100.000       ->  frame   2520 time 25200.000       ->  frame   2530 time 25300.000       ->  frame   2540 time 25400.000       ->  frame   2550 time 25500.000       ->  frame   2560 time 25600.000       ->  frame   2570 time 25700.000       ->  frame   2580 time 25800.000       ->  frame   2590 time 25900.000       ->  frame   2600 time 26000.000       ->  frame   2610 time 26100.000       ->  frame   2620 time 26200.000       ->  frame   2630 time 26300.000       ->  frame   2640 time 26400.000       ->  frame   2650 time 26500.000       ->  frame   2660 time 26600.000       ->  frame   2670 time 26700.000       ->  frame   2680 time 26800.000       ->  frame   2690 time 26900.000       ->  frame   2700 time 27000.000       ->  frame   2710 time 27100.000       ->  frame   2720 time 27200.000       ->  frame   2730 time 27300.000       ->  frame   2740 time 27400.000       ->  frame   2750 time 27500.000       ->  frame   2760 time 27600.000       ->  frame   2770 time 27700.000       ->  frame   2780 time 27800.000       ->  frame   2790 time 27900.000       ->  frame   2800 time 28000.000       ->  frame   2810 time 28100.000       ->  frame   2820 time 28200.000       ->  frame   2830 time 28300.000       ->  frame   2840 time 28400.000       ->  frame   2850 time 28500.000       ->  frame   2860 time 28600.000       ->  frame   2870 time 28700.000       ->  frame   2880 time 28800.000       ->  frame   2890 time 28900.000       ->  frame   2900 time 29000.000       ->  frame   2910 time 29100.000       ->  frame   2920 time 29200.000       ->  frame   2930 time 29300.000       ->  frame   2940 time 29400.000       ->  frame   2950 time 29500.000       ->  frame   2960 time 29600.000       ->  frame   2970 time 29700.000       ->  frame   2980 time 29800.000       ->  frame   2990 time 29900.000      Reading frame    3000 time 30000.000    ->  frame   3000 time 30000.000       ->  frame   3010 time 30100.000       ->  frame   3020 time 30200.000       ->  frame   3030 time 30300.000       ->  frame   3040 time 30400.000       ->  frame   3050 time 30500.000       ->  frame   3060 time 30600.000       ->  frame   3070 time 30700.000       ->  frame   3080 time 30800.000       ->  frame   3090 time 30900.000       ->  frame   3100 time 31000.000       ->  frame   3110 time 31100.000       ->  frame   3120 time 31200.000       ->  frame   3130 time 31300.000       ->  frame   3140 time 31400.000       ->  frame   3150 time 31500.000       ->  frame   3160 time 31600.000       ->  frame   3170 time 31700.000       ->  frame   3180 time 31800.000       ->  frame   3190 time 31900.000       ->  frame   3200 time 32000.000       ->  frame   3210 time 32100.000       ->  frame   3220 time 32200.000       ->  frame   3230 time 32300.000       ->  frame   3240 time 32400.000       ->  frame   3250 time 32500.000       ->  frame   3260 time 32600.000       ->  frame   3270 time 32700.000       ->  frame   3280 time 32800.000       ->  frame   3290 time 32900.000       ->  frame   3300 time 33000.000       ->  frame   3310 time 33100.000       ->  frame   3320 time 33200.000       ->  frame   3330 time 33300.000       ->  frame   3340 time 33400.000       ->  frame   3350 time 33500.000       ->  frame   3360 time 33600.000       ->  frame   3370 time 33700.000       ->  frame   3380 time 33800.000       ->  frame   3390 time 33900.000       ->  frame   3400 time 34000.000       ->  frame   3410 time 34100.000       ->  frame   3420 time 34200.000       ->  frame   3430 time 34300.000       ->  frame   3440 time 34400.000       ->  frame   3450 time 34500.000       ->  frame   3460 time 34600.000       ->  frame   3470 time 34700.000       ->  frame   3480 time 34800.000       ->  frame   3490 time 34900.000       ->  frame   3500 time 35000.000       ->  frame   3510 time 35100.000       ->  frame   3520 time 35200.000       ->  frame   3530 time 35300.000       ->  frame   3540 time 35400.000       ->  frame   3550 time 35500.000       ->  frame   3560 time 35600.000       ->  frame   3570 time 35700.000       ->  frame   3580 time 35800.000       ->  frame   3590 time 35900.000       ->  frame   3600 time 36000.000       ->  frame   3610 time 36100.000       ->  frame   3620 time 36200.000       ->  frame   3630 time 36300.000       ->  frame   3640 time 36400.000       ->  frame   3650 time 36500.000       ->  frame   3660 time 36600.000       ->  frame   3670 time 36700.000       ->  frame   3680 time 36800.000       ->  frame   3690 time 36900.000       ->  frame   3700 time 37000.000       ->  frame   3710 time 37100.000       ->  frame   3720 time 37200.000       ->  frame   3730 time 37300.000       ->  frame   3740 time 37400.000       ->  frame   3750 time 37500.000       ->  frame   3760 time 37600.000       ->  frame   3770 time 37700.000       ->  frame   3780 time 37800.000       ->  frame   3790 time 37900.000       ->  frame   3800 time 38000.000       ->  frame   3810 time 38100.000       ->  frame   3820 time 38200.000       ->  frame   3830 time 38300.000       ->  frame   3840 time 38400.000       ->  frame   3850 time 38500.000       ->  frame   3860 time 38600.000       ->  frame   3870 time 38700.000       ->  frame   3880 time 38800.000       ->  frame   3890 time 38900.000       ->  frame   3900 time 39000.000       ->  frame   3910 time 39100.000       ->  frame   3920 time 39200.000       ->  frame   3930 time 39300.000       ->  frame   3940 time 39400.000       ->  frame   3950 time 39500.000       ->  frame   3960 time 39600.000       ->  frame   3970 time 39700.000       ->  frame   3980 time 39800.000       ->  frame   3990 time 39900.000      Reading frame    4000 time 40000.000    ->  frame   4000 time 40000.000       ->  frame   4010 time 40100.000       ->  frame   4020 time 40200.000       ->  frame   4030 time 40300.000       ->  frame   4040 time 40400.000       ->  frame   4050 time 40500.000       ->  frame   4060 time 40600.000       ->  frame   4070 time 40700.000       ->  frame   4080 time 40800.000       ->  frame   4090 time 40900.000       ->  frame   4100 time 41000.000       ->  frame   4110 time 41100.000       ->  frame   4120 time 41200.000       ->  frame   4130 time 41300.000       ->  frame   4140 time 41400.000       ->  frame   4150 time 41500.000       ->  frame   4160 time 41600.000       ->  frame   4170 time 41700.000       ->  frame   4180 time 41800.000       ->  frame   4190 time 41900.000       ->  frame   4200 time 42000.000       ->  frame   4210 time 42100.000       ->  frame   4220 time 42200.000       ->  frame   4230 time 42300.000       ->  frame   4240 time 42400.000       ->  frame   4250 time 42500.000       ->  frame   4260 time 42600.000       ->  frame   4270 time 42700.000       ->  frame   4280 time 42800.000       ->  frame   4290 time 42900.000       ->  frame   4300 time 43000.000       ->  frame   4310 time 43100.000       ->  frame   4320 time 43200.000       ->  frame   4330 time 43300.000       ->  frame   4340 time 43400.000       ->  frame   4350 time 43500.000       ->  frame   4360 time 43600.000       ->  frame   4370 time 43700.000       ->  frame   4380 time 43800.000       ->  frame   4390 time 43900.000       ->  frame   4400 time 44000.000       ->  frame   4410 time 44100.000       ->  frame   4420 time 44200.000       ->  frame   4430 time 44300.000       ->  frame   4440 time 44400.000       ->  frame   4450 time 44500.000       ->  frame   4460 time 44600.000       ->  frame   4470 time 44700.000       ->  frame   4480 time 44800.000       ->  frame   4490 time 44900.000       ->  frame   4500 time 45000.000       ->  frame   4510 time 45100.000       ->  frame   4520 time 45200.000       ->  frame   4530 time 45300.000       ->  frame   4540 time 45400.000       ->  frame   4550 time 45500.000       ->  frame   4560 time 45600.000       ->  frame   4570 time 45700.000       ->  frame   4580 time 45800.000       ->  frame   4590 time 45900.000       ->  frame   4600 time 46000.000       ->  frame   4610 time 46100.000       ->  frame   4620 time 46200.000       ->  frame   4630 time 46300.000       ->  frame   4640 time 46400.000       ->  frame   4650 time 46500.000       ->  frame   4660 time 46600.000       ->  frame   4670 time 46700.000       ->  frame   4680 time 46800.000       ->  frame   4690 time 46900.000       ->  frame   4700 time 47000.000       ->  frame   4710 time 47100.000       ->  frame   4720 time 47200.000       ->  frame   4730 time 47300.000       ->  frame   4740 time 47400.000       ->  frame   4750 time 47500.000       ->  frame   4760 time 47600.000       ->  frame   4770 time 47700.000       ->  frame   4780 time 47800.000       ->  frame   4790 time 47900.000       ->  frame   4800 time 48000.000       ->  frame   4810 time 48100.000       ->  frame   4820 time 48200.000       ->  frame   4830 time 48300.000       ->  frame   4840 time 48400.000       ->  frame   4850 time 48500.000       ->  frame   4860 time 48600.000       ->  frame   4870 time 48700.000       ->  frame   4880 time 48800.000       ->  frame   4890 time 48900.000       ->  frame   4900 time 49000.000       ->  frame   4910 time 49100.000       ->  frame   4920 time 49200.000       ->  frame   4930 time 49300.000       ->  frame   4940 time 49400.000       ->  frame   4950 time 49500.000       ->  frame   4960 time 49600.000       ->  frame   4970 time 49700.000       ->  frame   4980 time 49800.000       ->  frame   4990 time 49900.000      Reading frame    5000 time 50000.000    ->  frame   5000 time 50000.000       ->  frame   5010 time 50100.000       ->  frame   5020 time 50200.000       ->  frame   5030 time 50300.000       ->  frame   5040 time 50400.000       ->  frame   5050 time 50500.000       ->  frame   5060 time 50600.000       ->  frame   5070 time 50700.000       ->  frame   5080 time 50800.000       ->  frame   5090 time 50900.000       ->  frame   5100 time 51000.000       ->  frame   5110 time 51100.000       ->  frame   5120 time 51200.000       ->  frame   5130 time 51300.000       ->  frame   5140 time 51400.000       ->  frame   5150 time 51500.000       ->  frame   5160 time 51600.000       ->  frame   5170 time 51700.000       ->  frame   5180 time 51800.000       ->  frame   5190 time 51900.000       ->  frame   5200 time 52000.000       ->  frame   5210 time 52100.000       ->  frame   5220 time 52200.000       ->  frame   5230 time 52300.000       ->  frame   5240 time 52400.000       ->  frame   5250 time 52500.000       ->  frame   5260 time 52600.000       ->  frame   5270 time 52700.000       ->  frame   5280 time 52800.000       ->  frame   5290 time 52900.000       ->  frame   5300 time 53000.000       ->  frame   5310 time 53100.000       ->  frame   5320 time 53200.000       ->  frame   5330 time 53300.000       ->  frame   5340 time 53400.000       ->  frame   5350 time 53500.000       ->  frame   5360 time 53600.000       ->  frame   5370 time 53700.000       ->  frame   5380 time 53800.000       ->  frame   5390 time 53900.000      
[DEBUG  ] 
[DEBUG  ] 
[DEBUG  ] GROMACS reminds you: "Drugs are Bad, mmokay" (South Park)
[DEBUG  ] 
[DEBUG  ] Note that major changes are planned in future for trjconv, to improve usability and utility.
[DEBUG  ] Select group for output
[DEBUG  ] Selected 20: 'GMXMMPBSA_REC_GMXMMPBSA_LIG'
[INFO   ] Building AMBER topologies from GROMACS files... Done.

[INFO   ] Loading and checking parameter files for compatibility...
[INFO   ] Preparing trajectories for simulation...

[INFO   ] 4 frames were processed by cpptraj for use in calculation.
[INFO   ] Running calculations on normal system...
[INFO   ] Beginning GB calculations with /home/romi/anaconda3/envs/gmxMMPBSA/bin/sander
[INFO   ]   calculating complex contribution...
[INFO   ]   calculating receptor contribution...
[INFO   ]   calculating ligand contribution...
[INFO   ] Beginning PB calculations with /home/romi/anaconda3/envs/gmxMMPBSA/bin/sander
[INFO   ]   calculating complex contribution...
[INFO   ]   calculating receptor contribution...
[INFO   ]   calculating ligand contribution...
[INFO   ] Parsing results to output files...

[DEBUG  ] Creating converter from 5 to 3
[INFO   ] Timing:
[INFO   ] Total GROMACS setup time:                   0.149 min.
[INFO   ] Total AMBER setup time:                     0.010 min.
[INFO   ] Creating trajectories with cpptraj:         0.003 min.
[INFO   ] Total calculation time:                     7.093 min.
[INFO   ] Total GB calculation time:                  1.035 min.
[INFO   ] Total PB calculation time:                  6.058 min.
[INFO   ] Statistics calculation & output writing:    0.003 min.
[INFO   ] Total time taken:                           7.260 min.

[INFO   ] 
   Finalizing gmx_MMPBSA: [ERROR  ] = 0; [WARNING] = 0
   Check the gmx_MMPBSA.log file for more details...

[INFO   ] 
 Thank you for using gmx_MMPBSA. Please consider supporting gmx_MMPBSA by citing our publication:
    Valdés-Tresanco, M.S., Valdés-Tresanco, M.E., Valiente, P.A. and Moreno E. 
    gmx_MMPBSA: A New Tool to Perform End-State Free Energy Calculations with GROMACS. 
    J Chem Theory Comput., 2021, 17 (10):6281-6291. Epub 2021 Sep 29. PMID: 34586825.
    https://pubs.acs.org/doi/10.1021/acs.jctc.1c00645

Also consider citing MMPBSA.py:
    Miller III, B. R., McGee Jr., T. D., Swails, J. M. Homeyer, N. Gohlke, H. and Roitberg, A. E.
    MMPBSA.py: An Efficient Program for End-State Free Energy Calculations.
    J. Chem. Theory Comput., 2012, 8 (9) pp 3314-3321

[INFO   ] Opening gmx_MMPBSA_ana to analyze results...

[INFO   ] Finalized...
